The Effects of the pep Nuclear Reaction and Other Improvements in the Nuclear Reaction Rate Library on Simulations of the Classical Nova Outburst

We have continued our studies of the Classical Nova outburst by evolving TNRs on 1.25Msun and 1.35Msun WDs (ONeMg composition) under conditions which produce mass ejection and a rapid increase in the emitted light, by examining the effects of changes in the nuclear reaction rates on both the observable features and the nucleosynthesis during the outburst. In order to improve our calculations over previous work, we have incorporated a modern nuclear reaction network into our hydrodynamic computer code. We find that the updates in the nuclear reaction rate libraries change the amount of ejected mass, peak luminosity, and the resulting nucleosynthesis. In addition, as a result of our improvements, we discovered that the pep reaction was not included in our previous studies of CN explosions. Although the energy production from this reaction is not important in the Sun, the densities in WD envelopes can exceed $10^4$ gm cm$^{-3}$ and the presence of this reaction increases the energy generation during the time that the p-p chain is operating. The effect of the increased energy generation is to reduce the evolution time to the peak of the TNR and, thereby, the accreted mass as compared to the evolutionary sequences done without this reaction included. As expected from our previous work, the reduction in accreted mass has important consequences on the characteristics of the resulting TNR and is discussed in this paper.Comment: Accepted to the Astrophysical Journa

Hydrogen-Accreting Carbon-Oxygen White Dwarfs of Low Mass: Thermal and Chemical Behavior of Burning Shells

Numerical experiments have been performed to investigate the thermal behavior of a cooled down white dwarf of initial mass M_{\rm WD} = 0.516 M_{\sun} which accretes hydrogen-rich matter with Z = 0.02 at the rate $\dot{M}=10^{-8}$ \msun \yrm1, typical for a recurrent hydrogen shell flash regime. The evolution of the main physical quantities of a model during a pulse cycle is examined in detail. From selected models in the mass range $M_{\rm WD} = 0.52\div 0.68$ \msunend, we derive the borders in the $M_{\rm WD}$ - $\dot{M}$ plane of the steady state accretion regime when hydrogen is burned at a constant rate as rapidly as it is accreted. The physical properties during a hydrogen shell flash in white dwarfs accreting hydrogen-rich matter with metallicities Z = 0.001 and Z = 0.0001 are also studied. For a fixed accretion rate, a decrease in the metallicity of the accreted matter leads to an increase in the thickness of the hydrogen-rich layer at outburst and a decrease in the hydrogen-burning shell efficiency. In the $M_{\rm WD}$-$\dot{M}$ plane, the borders of the steady state accretion band are critically dependent on the metallicity of the accreted matter: on decreasing the metallicity, the band is shifted to lower accretion rates and its width in $\dot{M}$ is reduced.Comment: 31 pages and 10 Postscript figures; Accepted for publication on Ap

Hydrodynamic models for novae with ejecta rich in oxygen, neon and magnesium

The characteristics of a new class of novae are identified and explained. This class consists of those objects that have been observed to eject material rich in oxygen, neon, magnesium, and aluminum at high velocities. We propose that for this class of novae the outburst is occurring not on a carbon-oxygen white dwarf but on an oxygen-neon-magnesium white dwarf which has evolved from a star which had a main sequence mass of approx. 8 solar masses to approx. 12 solar masses. An outburst was simulated by evolving 1.25 solar mass white dwarfs accreting hydrogen rich material at various rates. The effective enrichment of the envelope by ONeMg material from the core is simulated by enhancing oxygen in the accreted layers. The resulting evolutionary sequences can eject the entire accreted envelope plus core material at high velocities. They can also become super-Eddington at maximum bolometric luminosity. The expected frequency of such events (approx. 1/4) is in good agreement with the observed numbers of these novae

CNO abundances and hydrodynamic models of the Nova outbursts. 4: Comparison with observations

A variety of observations of novae are discussed in light of theoretical models. It is proposed that the nearly constant bolometric luminosity of FH Ser originates in the non-degenerate hydrogen-burning region at the bottom of the hydrogen-rich envelope which remains after the primary ejection. The shift of the wavelength of peak emission from the visual to shortward of the ultraviolet is caused by the decrease of the photospheric radius of the remnant envelope as the bolometric luminosity stays nearly constant. The oscillations in the light curve of GK Per during the transition stage can be explained by a pulsation of the remnant envelope when it is the size of the Roche lobe. The CNO over-abundances in novae reported by various observers are strongly suggestive of this nova mechanism. Finally, the implications of the upper limits of C-13 and N-15 in DQ Her are discussed

Numerical simulation of the magnetospheric gate model for X-ray bursters

A Lagrangian, fully implicit, one dimensional hydrodynamic computer code was used to investigate the evolution of a gas cloud impacting the surface of a 20 km, 1 Msub solar neutron star. This gas is initially at rest with respect to the surface of the neutron star, extends to 185 km above the surface, and is optically thick. The infall results in a burst which lasts about 0.1 seconds and reached a peak luminosity and effective temperature of 240,000 Lsub solar and 9 million; respectively. The burst was followed by a phase of oscillations with a period 0.2 seconds

A burst from a thermonuclear runaway on an ONeMg white dwarf

Studies which examine the consequences of accretion, at rates of 10(exp -9) solar mass/yr and 10(exp -10) solar mass/yr, onto an ONeMg white dwarf with a mass of 1.35 solar masses are performed. In these studies, a Lagrangian, hydrodynamic, one-dimensional computer code was used. The code now includes a network with 89 nuclei up to Ca-40, elemental diffusion, new opacities, and new equation of state. The initial abundance distribution corresponded to a mixture that was enriched to either 25, 50, or 75 percent in products of carbon burning. The remaining material in each case is assumed to have a solar composition. The evolution of the thermonuclear runaway in the 1.35 solar mass white dwarf, with M = 10(exp -9) solar mass, produced peak temperatures in the shell source exceeding 300 million degrees. The sequence produced significant amounts of Na-22 from proton captures onto Ne-20 and significant amounts of Al-26 from proton captures on Mg-24. This sequence ejected 5.2 x 10(exp -6) solar mass moving with speeds from approximately 100 km/s to 2300 km/s. When the mass accretion rate was decreased to 10(exp -10) solar mass, the resulting thermonuclear runaway produced a shock that moved through the outer envelope of the white dwarf and raised the surface luminosity to L greater than 10(exp 7) solar luminosity and the effective temperature to values exceeding 10(exp 7) K. The interaction of the material expanding from off of the white dwarf with the accretion disk should produce a burst of gamma-rays

Thermonuclear runaways in thick hydrogen rich envelopes of neutron stars

A Lagrangian, fully implicit, one dimensional hydrodynamic computer code was used to evolve thermonuclear runaways in the accreted hydrogen rich envelopes of 1.0 Msub solar neutron stars with radii of 10 km and 20 km. Simulations produce outbursts which last from about 750 seconds to about one week. Peak effective temeratures and luninosities were 26 million K and 80 thousand Lsub solar for the 10 km study and 5.3 millison and 600 Lsub solar for the 20 km study. Hydrodynamic expansion on the 10 km neutron star produced a precursor lasting about one ten thousandth seconds

Modeling CHANDRA Low Energy Transmission Grating Spectrometer Observations of Classical Novae with PHOENIX. I. V4743 Sagittarii

We use the PHOENIX code package to model the X-ray spectrum of Nova V4743 Sagittarii observed with the LETGS onboard the Chandra satellite on March 2003. Our atmosphere models are 1D spherical, expanding, line blanketed, and in full NLTE. To analyze nova atmospheres and related systems with an underlying nuclear burning envelope at X-ray wavelengths, it was necessary to update the code with new microphysics, as discussed in this paper. We demonstrate that the X-ray emission is dominated by thermal bremsstrahlung and that the hard X-rays are dominated by Fe and N absorption. The best fit to the observation is provided at a temperature of T_eff = 5.8 x 10^5 K, with L_bol = 50 000 L_sun. The models are calculated for solar abundances. It is shown that the models can be used to determine abundances in the nova ejecta.Comment: 8 pages, 6 figures, accepted for publication in Astronomy & Astrophysic